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Abstract 

A self-consistent mean-field theory explaining the thermal remagnetization (TR) of 
polycrystalline permanent magnets is given. The influence of the environment of a 
grain is treated by an inclusion approximation, relating the field inside the grain 
to the local field outside by means of an internal demagnetization factor n. For the 
switching fields and the fluctuations of the local fields around the mean field Gaus- 
sian distributions of widths o~ s and 07 resp. are assumed. The isothermal hysteresis 
curve, the recoil curves, and the TR in dependence on the model parameters n, a s , 
and o~f are calculated. Furthermore, the influence of the initial temperature and the 
strong dependence of the TR on the demagnetization factor of the sample are stud- 
ied, and it is shown that for reasonable parameter sets TR effects up to 100 % are 
possible. The theoretical results correspond well with the experimental situation. 
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1 Introduction 



It is common knowledge that increasing temperature destroys ferromagnetism, 
due to the related decrease of both the saturation magnetization and the co- 
ercivity. Thus an experiment, where an unmagnetized sample is heated to 
some hundred degrees and becomes thereby magnetized, is a challenge to our 
understanding of permanent magnetism. This effect was discovered in SmCos 
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about 25 years ago [1,2] and named thermal remagnetization (TR). The exper- 
imental procedure, sketched in Fig. 1, includes the saturation and isothermal 
dc-demagnetization via the points "1" and "2" and a subsequent heating. The 
largest TR-effect, i.e. more than 80 % of the saturation magnetization at the 
initial temperature, has been observed in well-aligned sintered SmCo 5 for a 
closed magnetic circuit [2,3]. In open circuits the TR depends strongly on 
the demagnetization factor N of the sample [4,5]. In well texturized NeFeB 
magnets the TR-effect is lower, but well measurable [5,6]. In contrast to the 
nucleation controled SmCo 5 - and NdFeB-magnets the effect is absent [4,7] or 
at least very small [8] in high coercivity Sm 2 Coi 7 sintered magnets, which are 
believed to be pinning controled. In spite of the good experimental situation 
the theoretical understanding is dissatisfying. There is agreement on the fact 
that the average coercivities of positively and negatively magnetized grains 
have to be different and that the temperature dependence of the coercivity is 
essential for understanding the TR. The latter was demonstrated by experi- 
ments on barium ferrite, where a sign-change of the temperature coefficient of 
the coercive field He results in an "inverse" TR effect, i.e. a remagnetization 
upon cooling [5]. Livingston and Martin [4] argued that the TR in SmCo 5 
is due to the change of the grains from single-domain to multi-domain state. 
Their mechanism can not explain TR to more than 50 %. We proposed an 
alternative mechanism [5], where the TR is caused by the fluctuating fields 
in the neighbourhood of "hard" grains. This first model aimed at a quali- 
tative understanding of the TR solely, therefore oversimplified distributions 
of the switching fields and the field fluctuations were used. In a subsequent 
calculation [9] we introduced Gaussian distributions both for the switching 
fields and for the field fluctuations. Furthermore the distribution width of the 
latter was chosen as function of the average magnetization. This brought the 
resulting TR curves closer to the experimental situation, especially well below 
the temperature T max , where the maximum TR occurs, but was not able to 
explain the vanishing TR at temperatures above. The same kind of distribu- 
tions for the fluctuation fields were used by Miiller et. al. [6]. They presented 
a mean- field theory devoted mainly to the coercivity of NdFeB, where non- 
magnetic phases as well as multi-domain magnetic grains are present. They 
also derived a formula for the maximum TR. Whereas they took the field 
fluctuations into account while calculating the demagnetization curve, they 
neglected these fluctuations in that part of their paper, which was devoted to 
the TR. The resulting overestimation of the TR for NeFeB they attributed 
to that neglection of the field fluctuations, and furthermore to the pinning of 
domain walls, as well as to the influence of the sample demagnetization factor. 
Supported by the experimental fact, that SmCo 5 of lower to medium coerciv- 
ity in its dc-demagnetized state contains multi-domain grains [10,4,11], Lileev 
et al. [12,13] presented a theory, which takes into account both multi-domain 
grains, and field fluctuations caused by the interacting field of neighbouring 
grains. Since they included the mutual dipol-dipol interaction of the grains 
they had to carry out numerical simulations which makes it difficult to handle 
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their theory, especially if one is interested in the influence of basic parameters 
e.g. the distribution width of the field fluctuations, the influence of the ini- 
tial temperature etc.. They calculated a TR curve which has a rather sharp 
peak. Whereas the magnitude of the maximum effect is comparable to the 
experiment, the shape of the theoretical curve differs from the experimental 
situation. Nevertheless, we share the opinion that the influence of the nearest 
neighbourhood of a grain is an essential for every elaborated theory of the 
TR. For this task the current paper extends our former theory [9] by an in- 
clusion approximation. Furthermore we take into account the possibility that 
the grains may nucleate into multi-domain states. The paper is organized as 
follows: In the next section we present the model. Since a theory of the TR 
presumes a calculation of both the isothermal hysteresis loops and the recoil 
curves from every point of the hysteresis loop, we present the related theory 
and the results in subsection 3.1 and 3.2 respectively. The calculation of the 
TR due to heating is the topic of subsection 3.3. In the last section we discuss 
the results. 



2 The model 

The polycrystalline permanent magnet is an ensemble of high-coercive mag- 
netic uniaxial grains. To fix our model we have to define the properties of both 
a single grain and the ensemble. 

2. 1 Properties of a single grain 

Every grain is characterized by 

(Al) Ms(T), the temperature dependent saturation magnetization, 

(A2) H S (T), the temperature dependent switching field. This is the absolute value 

of the field needed for reversing the magnetization of the grain in a closed 

circuit. 

(A3) n, the "internal" demagnetization factor, which is determined by the shape 
of the grain, 

(A4) c, the direction of the easy axis due to the high uniaxial anisotropy 
(A5) V, the volume of the grain. 

2.2 Properties of the ensemble 
Regarding the ensemble we assume 
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(Bl) all grains exhibit the same Ms(T), 

(B2) the easy axis c of the grains are completely aligned (ideal texture), 

(B3) the grains differ in their switching fields H s , but the distribution of these 

switching fields should be known. For simplicity we assume a Gaussian 

normalized with respect to the region < H s < oo, 



with the "mean switching field" H s and the distribution width a s . 
(B4) the temperature dependence (not the absolute value!) of the switching fields 
of all grains is equal. 

It is problematic to get direct information both on the temperature depen- 
dence of the switching fields and their distribution. What can be measured 
is the average magnetization in dependence on the external magnetic field, 
the temperature and the time. The switching field distribution can not be 
measured directly. In non-interacting ensembles the switching field distribu- 
tion may be determined from remanence measurements [14] , but if the grains 
interact strongly this method fails. From the Henkel plots measured for Nd- 
FeB [15,16] we know that most sintered hard magnets are strongly interacting 
ensembles. Thus, for the moment, we have to regard H s and o~ s as model pa- 
rameters. Fortunately we will see later on that our theory delivers a possibility 
to relate H s to the coercivity H c of the sample. 

2.3 Inclusion approximation 

Since we believe that the TR is caused mainly due to grain interactions we 
want to explain the approximations in detail. On a macroscopic length scale 
the internal magnetic field and the magnetization of the sample are homoge- 
neous. But, if the length scale is reduced to the order of some grain diameters 
the magnetization becomes coarser and the inhomogeneities gain influence. 
Averaging over volume elements containing a few grains only, yields values, 
which deviate from the average. Thus, we consider every grain as an inclusion 
embedded in a local environment, which may differ stochastically in its mag- 
netic field and magnetization from the related averaged values. We show this 
schematically in Fig. 2. For simplicity we approximate the magnetization in the 
vicinity of the inclusion by the mean magnetization (M), neglecting the men- 
tioned fluctuations. Otherwise we allow for field fluctuations AH = H — (H) 
around the mean internal field in the environment of a grain. The mean inter- 
nal field is related to the external applied field H ext according to 




(1) 



fi (H)=fi H ext -N{M) 



(2) 
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with N being the demagnetization factor of the sample. The field fluctuations 
are characterized by the following assumptions 



(CI) The local magnetic fields (in c-direction) are Gaussian distributed 
1 -(^t) 2 

f(H) = c V */ J with AH = H-(H) . (3) 

(C2) The distribution width &f is itself a function of the mean magnetization. 
We set 




*f = ff f\ 1 - 1 7Pr) ■ ( 4 ) 



Assumption CI is a consequence of the central limit theorem of probability 
theory. Of course the local fields fluctuate around the mean internal field also in 
direction, thus, strictly speaking eq. (3) accounts for the fluctuations of the z- 
component only. Furthermore the deviation of the field direction from the easy 
axis gives rise to rotation processes, but this is surely a small effect in ideally 
aligned magnets with high anisotropy constants. Assumption C2 attempts to 
take into account that in an ideally textured magnet the fluctuation width 
is zero in the saturated state and maximum in the dc-demagnetized state. If 
cff((M)) is expanded to second order and the expansion parameters are fixed 
by these requirements together with the demand that <7/((M)) has to be a 
symmetric function one gets eq. (4). 



3 Field- and temperature-dependent magnetization 



3. 1 Calculation of the isothermal demagnetization curve 



After the saturation all grains are up (f) magnetized. They posses their max- 
imum switching fields and (M) = M s holds, due to assumption B2. If the 
external field is lowered to H^J t (point "1" in Fig. 1) the magnetization of a 
grain is switched down (J,) if its internal field H$ is lower than — H s and the 

internal field after switching will be smaller than +H S . Thus we have the 
two conditions 

H$<-H a ("switching condition") , (5) 
H^<+H s ( "not-back-switching condition" ) . (6) 

The internal fields before and after switching are with respect to the above 
introduced approximations: 
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VoHij =^oH — n(M s — (M)i) , (7) 
f i H^= f i H-n(-M s -(M) 1 ) . (8) 

Here H is the local field in the environment of the grain. A variation of the 
external field H ext results in a change of (M) on the one hand and on the other 
hand it changes the probability f(H)dH that the field in the environment of 
a grain is between H and H + dH due to the relation (2) and eq. (3). If the 
switching condition (5) is fulfilled but the condition (6) is not fulfilled, i.e. 



H$>H a ("back-switch condition") (9) 

holds, the grain cannot jump into a stable state. If the grain is large enough it 
solves the conflict by an incomplete jump, what turns its single-domain state 
(SDS) into a multi-domain-state (MDS). We will call such grains "weak". 
Whereas the "hard" grains can exist in states with ±Mg only, the magneti- 
zation of a weak grain (M)j is an average over the upwards and downwards 
magnetized volume fractions within the grain. From phase theory of 180°- 
domains we get for the averaged magnetization of the i th grain in dependence 
on the local field in the environment (cf. Appendix A) 



(M)i = (M) + l^ . (10) 
n 

Opposite to the hard grains the weak grains have no memory since the mag- 
netization follows the local field immediately as long as |(M)j| < Ms holds. 
Otherwise they are saturated up- or downwards. If one calculates the proba- 
bility that a grain with a given H s is in a MDS the switching condition (5) 
limits the integration over the local field distribution from above, whereas the 
back-switching condition (9) limits from below. The related probability is 



Pw (H s )= JdHf(H) (11) 
The two limiting fields are 



fi H H = -fi H s + nM s - n(M) , (12) 
fi H L = u, H s - nM s ~ n(M) . (13) 

The requirement that the upper limit has to be greater then the lower one 
yields the condition 

fi H s < nM s . (14) 
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Thus, grains with switching fields smaller than their own internal demagne- 
tizing field are weak. The magnetization of the volume fraction of the weak 
grains with a given H s is 



(M) W (H S ) = M S J dHi{H)-M s J dHi(H)+ J dHi(H)(M) i .(15) 



H H -oo H L 

The integral may be evaluated yielding 



(M) W (H S ) = (erf(y H ) + erf(y L )) (16) 



+ \ ( ^ + (M)i(l ~ ] (er%H) - erf(y L )) 



2 \ n ^ n 



with 



(e-v 2 H - e -vVj 



-fi H s + nM s - n{M) l - fi Hjl\ + N(M) 1 
VH Wf (1 - (M)j/Ml) [ n 



fi H s - nM s - n(M) 1 - fipH^ + N(M) l 
1 W (1-(M)?/M|) 



The index 1 at H^J t and (M)i indicates that the calculation is for point 1 in 
Fig. 1. 



For hard grains is higher than Hh- If these grains fulfill the switching 
condition (5), condition (6) will be fulfilled automatically. The probability to 
find a grain downwards magnetized if a field H^) t is applied is 



-t-cxj 

Pl (H s ) = JdH f (1) (tf) &(-H s - H$) &(H S - H$) , (19) 

— oo 

with and H-^ being the fields within the grain in the upwards and down- 
wards magnetized state resp. and Q(x) is the Heaviside function. f(i) is the 
field distribution eq. (3) with Hq{H) = fipH^Jt — N(M)i. The integration in 
eq. (19) is easily done yielding 



Pl (H s ) = ^(l + eri(x H )) (20) 
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. + , -VoH s + nM s - n{M) l - ^ggg + N(M) 1 

With Xh = / ,,i9 / ,, 9N 

(MO<T f (1 - (M)?/M|) 

For the magnetization of such a fraction of hard grains with a given H s we 
find 



(M) 1 (H S ) = M S (1 - 2 Pl (tf s )) = -M 5 erf(a^) (21) 

The total magnetization of the sample one gets from averaging with respect 
to 



nMs/no oo 

(M>!= | dH s g(H s )(M) wA (H s )+ J d H s g(H s ) (M)i(H s ) (22) 

nM s /no 

Here the integration has been splitted due to the different contributions of the 
weak and hard grains. For the weak grains we have to use eq. (17) with H ext 
being H^ x \ and (M) being (M) 1 . Since the right hand side of eq. (22) depends 
on the mean magnetization {M)\ itself, we have to solve this implicit equation 
numerically. This is done by fixing (M)i/Ms and searching the related H^ x \. 
Fortunately there is always exactly one solution for —1 < (M)i/Ms < 1. 
The aim of the present paper is to clarify, how the TR depends on the model 
parameters a s , Of as well as on the internal demagnetization factor n. We vary 
these parameters around the values M s = IT, a s = 1.5T//x , oj = 0.5T//i , 
n = 0.333, and a mean value of the switching field distribution of ixqH s = 3.0T. 
The values of Ms, H s , and a s resemble SmCos (VACOMAX 200), where we 
measured at T=300K e.g. Ms = 0.972T, ^qHq = 2.87T. a s is a rough estimate 
of the width of the differentiated experimental demagnetization curve. The 
results for the calculated demagnetization curve are plotted in Fig. 3, where 
a spherical sample was assumed (N = 1/3). 



3.2 Calculation of the isothermal recoil curve 



Next, as shown in Fig. 1 (Point "1" to point "2"), the external field is changed 
from H^Jt to H^ t with H^ t > H^) t via the recoil curve. Some of the hard 
"down" -grains are switched back, whereas the weak grains shift their magne- 
tization reversible. The related mean magnetization is (M) 2 . The contribution 
of the weak grains results from eq. (17), if we insert HfJ t for H ext and (M} 2 
for (M) resp. followed by an integration from zero to nMs/ Ho with respect to 
H s . The contribution of the hard grains is more difficult to calculate, due to 
the memory effect. Let us consider the probability w^{H s ) that a hard grain 
with a given H s was switched by H^ x \ and switched back by Hf x \ afterwards. 
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To calculate this probability it is inevitable to make an assumption on the 
correlation between the fluctuation of the local field Hi which acts on a grain 
if H^rf is applied and the fluctuation of H 2 according to the external field 
H^} t . Of course, if H^t is only slightly different from H^ t it is unlikely that 
the neighbourhood of a grain changes considerably, so that Hi and H 2 should 
be strongly correlated. With increasing distance between H^) t and H^ t this 
correlation will vanish, due to the multitude of switching processes. The rea- 
soning is as follows: The thermodynamical potential has surely a lot of nearly 
degenerated local minima. Applying H^ t selects one of them. A small change 
of the applied field is not enough to overcome the barrier between adjacent 
minima, but larger changes generate energy surfaces, which are completely 
different in their topological structure, thus making states accessible, which 
may be very different from the state at H^) t . The change to such a state will 
be accompanied by a lot of correlated switchings. Even if the mean magneti- 
zation is changed only slightly, the neighbourhood of an individual grain may 
be completely different, hence we can average with respect to Hi and H 2 in- 
dependently, if the difference H^) t — H^ t is not to small. Thus the probability 
Wtf decouples accordingly: 

w u (H s )=pi(H s )q 2 (H s ) (23) 
withpi(if s ) from eq. (21) and 

oo 

q 2 (H s )= J dHi {2) (H)e(H%-H s )e(H% + Hs ) 

— oo 

= \ (l - erf(y L )) (24) 

with 



= fM>H a ~ nM s - n(M) 2 - fi H ( e 2 J t + N(M) 2 
VL WTf{l-(M)l/M?) 

Again the switching conditions are included by help of related Heaviside func- 
tions. Due to the mentioned magnetization changes in the environment of a 
grain it may happen that grains which resisted H^) t are switched down by 
HfJt ("up"- "down" -contribution). The related probability factorizes also 

w u (H s ) = (l-pi(H s ))p 2 (H s ) (26) 
with pi(H s ) again from eq. (21) and 
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p 2 (H s ) = J dHf {2) (H)Q(-H s - H®) Q(H S - H%) (27) 

— oo 

= \ (l + ed(z B )) (28) 



with 



-fi H s + nM s - n(M) 2 - /xpgg + N(M) 2 
z h = T, 777727777S ■ ( 29 ) 



fJ, (Tf 



(l - (M) 2 2 /M; 



There are two further probabilities corresponding to the remaining two histo- 
ries^ which a grain may experience, i.e. the probability that it resisted both 
H^xt and H^ x \ ( "up" -"up" contribution) and the probability that a grain was 
switched by H^ x \ but resisted back-switching by Hf x \ ( "down" - "down" contri- 
bution). We find 

w u (H s ) = (l- Pl (H s ))(l-p 2 (H s )) , (30) 
w ll (H 9 )=p 1 (H a )(l-q 2 (H a )) . (31) 

Thus hard grains with a given H s contribute to the magnetization 

M 2 (H s ) = M s (w n (H s )-w n (H s )+w^(H s )-w n (H s )) . (32) 

Averaging with respect to H s yields the recoil curve for the mean magnetiza- 
tion (M) 2 in state 2 in dependence on both magnetic fields H^ x \ and H.f x \: 



nM s /n oo 

(M) 2 = j dH s g(H s )(M) w , 2 (H s )+ J dH s g(H s )(M) 2 (H s ) (33) 

nM s /no 

In Fig. 3 we show additional to the recoil curve a set of minor loops starting 
from different values of (M) 1 . 



3. 3 Calculation of the thermal remagnetization 



In the preceding sections we calculated the demagnetization and the recoil 
curves with the implicit understanding that the magnetization changes isother- 
mally at a temperature T . The actual TR occurs if a dc-demagnetized RE- 
TM-magnet is heated while Hf x \ is kept constant (cf. Fig. 1 from point "2" to 
point "T"). Raising the temperature has usually two effects. On the one hand 
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the saturation magnetization decreases with increasing temperature and on 
the other hand the switching fields H s are changed. Whereas the temperature 
dependence of the saturation magnetization can be measured easily, it is im- 
possible to measure it for the switching fields. What can be measured directly 
is the temperature dependence of the coercive field. In appendix B we show 
that in our model H s is related to He and Ms by the simple formula 



if the influence of the weak grains is negligible. Otherwise we have to solve 
eq. (22) for H s with H^J t = H c and (M) 1 = 0. H C (T) and M S {T) are taken 
from measurements [17]. Due to the lack of better information we adopt the 
temperature dependence of H s for arbitrary values of the switching fields 



The calculation of the magnetization (M) 2 at an enhanced temperature T 
is similar to the calculation of the magnetization (M) 2 by help of eq. (33), 
whereby the deformation of the switching field distribution has to be regarded. 
That should be done by indicating the changed parameters with an index T, 
whereas the related parameters at the temperature To will not be marked. 
From eq. (35) and the normalization of the switching field distribution follows 
that a s has the same temperature dependence like H S (T). The number of weak 
grains increases with temperature, since H S {T) is usually much more reduced 
than nMs(T) if the temperature is raised. The magnetization of a weak grain 
with given H s at the temperature T results from eq. (17) with the appropriate 
H S (T), M S (T), and (M} 2 inserted. For the hard -grain contribution we have 
to calculate the probabilities according to eqs. (26,23,30,31), however for the 
elevated temperature. We find for the magnetization of the hard-grain fraction 
with given H s 



M^(H S ) = M T s (w^ T (H s ) - w n T (H s ) + w u T (H s ) - w u T (H s )) (36) 



It is obvious that the probability pi(H s ) remains unchanged, since the grains 
were switched downwards at the initial temperature T . Therefore we have 



fi H s (T)= f i H c (T) + nM s (T) 



(34) 



H S (T) = H S (T ) 



H S {T) 
H S (T ) 



(35) 




(37) 
(38) 
(39) 
(40) 



w u T (H s )= Pl (H s )qT(H s ) 
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The probabilities 



\ (l - erf (2/ 



and 



(41) 



f£(#.) = \ (l + erf(4)) (42) 

become temperature dependent both directly due to H S (T), M S (T), and 07 (T) 
and indirectly via (M) 2 , since we have now 



Vl — 



n(M) T 2 - pL Hil\ + N(M) 



(43) 



and 



z t = -VoHl + nM T s - n{M) T 2 - ^gg + N(M) T 2 
Hoaj(l-((Mf 2 /Mjy) 

The temperature dependence of aj = <7f(T) was assumed to be that of M s , 
since the fluctuations are caused by the inhomogeneities of the magnetization. 
From the numerical solution of the implicit self-consistent equation we calcu- 
lated the TR in dependence on the parameters for the above mentioned pa- 
rameter set. Here the starting temperature was fixed to 300 K and we assumed 
a spherical sample. The influence of the switching field distribution width is 
studied in Fig. 4 A. To get noticeable TR-effects the distribution width a s 
should not be too small. The dependence of the TR on the fluctuation width 
a/ is shown in Fig. 4 B. The TR peak becomes broader and shifts to lower 
temperatures with increasing field fluctuations. The strong dependence of the 
TR on the internal demagnetization factor is depicted in Fig. 4 C. For small 
demagnetization factors, i.e. for elongated grains, the TR is small, whereas 
effects up to 100 % are possible for platelet-shaped grains. Up to now there 
are no experimental results available regarding the influence of the internal 
demagnetization factor n onto the TR. Regarding the external demagnetiza- 
tion factor iV it is known that in closed-circuit measurements TR-effects up 
to 100% are possible [2] (related to M s (T max ) !) in the case of SmCo 5 sintered 
magnets. For samples measured in a VSM, where the demagnetization factor 
is between 0.1 and 0.6, a significant lowering of the TR with increasing demag- 
netization factor N is reported. This is exactly what the calculated curves in 
Fig. 5 demonstrate. Another feature, which has been studied experimentally, 
is the dependence of the TR on the initial temperature T . The maximum 
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TR increases and the peak position shifts slightly to lower temperatures with 
decreasing T [17], what agrees well with our theoretical results depicted in 
Fig. 6. 



4 Discussion 

The above presented theory explains the observed features of the TR-experi- 
ments well, whereby we focused especially to the shift of the temperature T max , 
where the maximum TR occurs, as well as to the height of the maximum it- 
self. Especially the strong dependence on the external demagnetization factor 
and on the starting temperature can be explained. Furthermore the theory is 
able to explain TR effects of more than 50 %. To proceed to the description 
of the TR it was inevitable to develop a theory for the hysteresis and the 
recoil curves. Whereas there exist a multitude of micromagnetic calculations 
trying to explain the coercivity by assuming special microstructures which 
give rise to either pinning of domain walls or nucleation, our theory starts a 
little above this level, since we do not judge about the reason for a switching 
field, but simply accept its existence. In this point our theory resembles the 
Preisach [18] model, which was created to take into account the interaction 
between different grains (bistable magnetic units). The main disadvantage of 
the Preisach model is due to the fact that the Preisach distribution function 
is not dependent on the magnetization state [18,19]. It is easy to prove that 
a theory of the TR based on this model is not able to explain TR effects of 
more than 50%. The reason is the lack of feedback. Nevertheless, we share the 
opinion that the difference between the local fields and the mean field, which 
is due to grain interactions, is the crucial point. In this work the width of these 
fluctuations is given by <j/. Our results demonstrate that without local field 
fluctuations there will be no TR effect above 50 %, e.g. in a normal mean-field 
theory. To take into account the feedback of the neighbourhood of a grain the 
inclusion approximation, albeit at a very rough level, was introduced. This 
provides us with a new parameter n, i.e. the internal demagnetization factor. 
Assigning the same internal demagnetization factor n to all grains seems to 
be a severe simplification, but since we have no reliable information on the 
switching field distribution, we probably account, at least partly, for different 
grain shapes and/or exchange coupling via the grain boundaries by adjusting 
the switching field distribution. Thus, n is more a feedback parameter than a 
real demagnetization factor. The most difficult problem is the determination 
of the temperature dependence of the switching fields. All the former theories 
used the temperature dependence of the coercivity for this task, due to the 
lack of better information. But it is evident that this forces the TR magnitude 
to drop to zero at the temperature Tg c=0 where the coercivity vanishes. Con- 
trary, investigations on SmCo 5 demonstrated that the TR manifestly remains 
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above Th c =o [17]. This is due to the fact that for this compound the difference 
between T# c=0 and the Curie temperature Tq is about 250 K [20]. Since our 
theory relates the mean switching field to the coercivity, it was possible to 
extract the temperature dependence H S (T) from the measured H C (T) curves. 
The derived formula is in congruence with the empirical formula used in [21,22] 
if one relates cHa (c^C 1) to the switching field and n to the corresponding ef- 
fective demagnetization factor N e ff. This holds as long as the mean switching 
field is greater than the internal demagnetizing field. Otherwise the simple rule 
given in eq. (34) fails. Another interesting point is whether the weak grains are 
essential for the description of the TR. The calculated relative volume fraction 
of weak grains in dependence on the temperature is depicted in a dashed fash- 
ion in Fig. 8. At T max the H s equals roughly nMs/i^o, i-e. the maximum TR is 
observed shortly before the majority of the grains become weak. At tempera- 
tures higher than T max the magnetization decreases. Sm-Co-magnets undergo 
irreversible changes in their microstructure at higher temperatures. That is 
why measurements above T max are rather difficult and, therefore, seldom re- 
ported. For the parameter set used in Fig. 7 the influence of weak grains is 
small as long as the temperature is below T max . For T > T max the neglection of 
the weak grains would result in a qualitative different behaviour, since we get 
a sign change of the TR before it finally vanishes, as shown in Fig. 7 in dotted 
manner. The reason for that peculiarity is made clear in Fig. 8, where we have 
plotted the temperature dependence of the probabilities w^, w^, w^, and 
Wn regarding the hard grain fraction. It is obvious that at first increases 
and when it drops down rises up, thus explaining the sign change. This 
behaviour is not observed if the probability p 2 , i.e. that a grain which resisted 
H^ x \ will be switched by H^ x \, is neglected, as it was done in former theories, 
but then an overestimation of the TR occurs, as visible in Fig. 7. For the de- 
scription of the TR at higher temperature it is necessary to take into account, 
that above T max even the switching fields of the hardest grains are lowered 
enough to be switched by the still existing stray fields. The latter do not van- 
ish until the magnetization breaks down. It is worth to point out, that weak 
grains are not necessary to get effects of more than 50 %. Such a large TR one 
gets due to the feedback of the magnetization via the internal demagnetizing 
fields. In some sense this models the well known avalanche effects observed 
during hysteresis measurement, which should also occur to a certain extent if 
the sample remagnetizes. In consideration of the various approximations one 
can not expect that all magnetic measurements are explained quantitatively 
by the presented theory. Nevertheless in Fig. 9, where we present a fit of the 
theory to our measurement of the TR in a SmCos sample (VACOMAX 200), 
the agreement is rather well, especially in comparison with Fig. 2 in ref. [13]. 
If one calculates the hysteresis curve with the parameters fitted to the TR 
curve the agreement is not as well as for the TR. This is not astonishing, 
since the theory was designed mainly for the qualitative understanding of the 
TR, which is very sensitive to the field fluctuations but less sensitive to the 
switching field distribution. For the hysteresis the situation is vice versa. To 
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get a good description of the hysteresis the H s spectrum has to be determined 
more accurately. Another reason may be the simplification in the description 
of the grain interaction. Correlated switching of a multitude of grains is taken 
into account in a simplified manner via n only, as discussed above. Surely 
this is more important for the hysteresis than for the TR. The neglection of 
the angular fluctuations of the field and the reversible rotation processes are 
in this regard of less importance, since this will influence both the hysteresis 
and the TR similarly. Since the present paper is focused on the theory, we 
save the detailed analysis of experimental results, as well as the modifications 
necessary to describe multi-phase magnets and the inverse TR in Ba-ferrites 
to forthcoming papers. In conclusion we may say that the presented theory 
is able to describe the rather complex findings by help of only three internal 
parameters, i.e. a s , <jf, and n, which may be helpful in characterizing different 
hard magnetic materials more effectively. 
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Appendices 

A The magnetization of a weak grain 

We want to find the mean internal magnetization of a grain in a MDS. For 
that purpose we assume, that the grain contains a lot of domains separated 
by 180°-domain walls. The dimension in field direction is large with respect 
to the dimension in a perpendicular direction. Thus the demagnetizing factor 
of such a slap-shaped domain is approximately zero. The average internal 
magnetization is then 



with A being the relative volume of the downwards magnetized domains. The 
magnetic energy density contains both the Zeeman energy 





w H = -{M) i H 
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and the energy of the stray fields 



IV stray 



J H dem d(M) t 



(A.3) 



Substituting for the demagnetization field 



n((M) i -(M)) 



(A.4) 



an integration of eq. (A.3) yields the total energy density w to t in dependence 
on A to be 



If the A, which minimizes the energy, is inserted into eq. (A.l) one gets eq. 
(10), as long as A is between zero and one. Otherwise the grain is saturated in 
forward or backward direction. Thus the MDS grains show no memory effect, 
but follow the local field immediately. 



B Relation between H c and H s 

The influence of the weak grains is negligible, if fioH s ^> nMs holds and the 
distribution width a s is small enough. The first integral in eq. (22) may be 
neglected and in the second one we can extend the lower integration limit 
to minus infinity. If the external field H^) t equals —H c we have (M)i = 
simultanously. Therefore we find now 



Substituting x = (H s — H s )/a s the switching field distribution becomes sym- 
metrical with respect to x — 0, thus the monotonous function erf(xfr) has to 
be asymmetrical to reduce the integral to zero. Otherwise we find from eq. (21) 
that the erf-function is shifted, since we have xh = {—H s — a s x + nMs/ '//o + 
Hc)/<Jf. To make the erf-function asymmetrical, the shift has to be zero. This 
yields eq. (34). 




(A.5) 




(B.l) 



— oo 
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Captions of figures 



Fig. 1. Scheme of the TR-experiment. From the saturated state the sample is 
isothermally driven to point "1" and afterwards along the recoil curve to point "2" . 
The TR happens, if the sample is heated (point "2" to point "T"), while the external 
field is kept constant (enhanced insert on the right hand side) . For the external fields 
at point "1" and "2" the sheared remanence coercivity H^ 1 and H ext = zero resp. 
were chosen, since the most TR-experiments start from the dc-demagnetized state. 

Fig. 2. The grain is considered as an inclusion with magnetization ±M$ (here we 
assumed —Ms, symbolized by the long white arrow) embedded in an environment 
differing in its magnetic field H from the mean magnetic field (H) within the sample 
due to the fluctuation AH. Fluctuations of the magnetization in the environment 
are neglected, thus the magnetization in the environment is equal to the mean 
magnetization of the sample. This is symbolized by the two shorter white arrows of 
equal length. 

Fig. 3. The demagnetization curve, the recoil curve, and a set of minor loops starting 
from different values of (M) 1 calculated for a representative parameter set close to 
SmCos-mangets at 300K. Both the sample and the grains are adopted as spheres. 

Fig. 4. The TR in dependence on (A) the width of the switching field distribution 
a s , (B) on the width of field fluctuations aj, and (C) on the internal demagnetization 
factor n. 

Fig. 5. The TR for different external demagnetization factors N. For simplicity this 
figure was calculated with H S {T) = Hc{T). 

Fig. 6. The TR for different values of the initial temperature To (parameters on 
the curves). 

Fig. 7. The TR for different levels of approximation. The solid line gives the results 
of the presented theory. The dot-dashed line results from neglecting the weak grain 
fraction and for the calculation of the dashed line, both the weak grains and p% , i.e. 
the probability, that a grain which resisted H^) t will be switched while heating at 

(2) 

Hext = 0' were neglected. 

Fig. 8. The probabilities w^ T , w ^ T , w^ T , and wn T (cf. eqs. (37-40)) in dependence 
on the temperature. The dashed line shows the ratio of the volume of weak grains 
to the total volume. 

Fig. 9. The TR and H C (T) measured for a SmCo 5 sample (VACOMAX 200) to- 
gether with the calculated TR, where the parameters a s , aj, and n have been 
adjusted accordingly. 
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